function model=update_zeta_2(model)



temp      = (sum(exp(model.alpha_2 + (1/2)*model.sigma_2.^2), 2))';
ind       = find(abs(temp)<model.MINVALUE);
temp(ind) = model.MINVALUE;

model.zeta_2 = temp;


end


